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ABSTRACT: Heat transfer coefficients of dryers are useful tools for correlation formulation and performance 
evaluation of process design of dryers as well as derivation of analytical model for predicting drying rates. A 
model equation for predicting heat transfer coefficient of air in a batch dryer using Buckingham Pi-theorem and 
dimensional analysis at various velocities has been formulated. The model was validated by drying unripe 
plantain chips in a batch dryer at air velocities of0.66, 0.92 and 1.20m/s at temperatures of42, 54.45 and 66°C, 
respectively. Based on the analogy of heat and mass transfer rate equations for constant drying period, the 
prediction from the developed model reasonably agreed with the experimental data. 
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I. INTRODUCTION 

Drying is a kinetic process that involves the removal of liquid, usually water from a moist material: solid, liquid 
or gas. The use of heat to remove liquid distinguishes drying from mechanical methods of removing water, such 
as: centrifugation, decantation, sedimentation and filtration in which no change in phase from liquid to vapour is 
experienced. The application of heat to remove moisture can also be found in processes like distillation and 
absorption. In the drying of solids to remove water a specialized device called dryer is used, and the desirable end 
products are in solid form. The final moisture contents of the dried solids are usually less than 1%. The chemistry 
of drying a moist material can be represented as: 

Moist material _solid + vapor (1.1) 

When the heat transfer by pure convection is used to dry a wet solid, the heat supplied is solely by sensible heat 
in the drying gas stream. A dynamic equilibrium exist between the rate of heat transfer to the material and the rate 
of vapor (mass) removal from the surface at instance, (that is, drying rate) and may be represented as follows. 


dx . AT 

— - hA - 

dt A 


( 1 . 2 ) 


The area of the heat and mass transfer may be assumed to be approximately equal [1], The study of convective 
heat transfer is centered on ways and means of determining the heat transfer coefficient, h for various flow regions 
(laminar, transition or turbulent flow) and over various geometries and configurations. The local and average heat 
transfer coefficient may be correlated by equations (1.3) and (1.4) respectively 

Nu x = f (x*, Re,, P r ) (1.3) 

Nu x = f (Re,, P r ) (1.4) 


where the subscript x emphasize the condition at a particular location on the surface. The problem of convection 
involves how these functions are obtained, there are two approaches: theoretical and xperimental. Theoretical 
approach involves solving the boundary layer equation for a particular geometry and equation such as equation 
(1.5) 


Nu x 


hL ST * , 

— = + - 1 , 

k ~ 8y* y 


(1.5) 
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which is a dimensionless temperature gradient at the surface. In the experimental approach, for a prescribed 
geometry in a parallel flow, if heated, convection heat transfer coefficient which is an average associated with the 
entire system could then be computed from Newton’s law of cooling. And from the knowledge of the characteristic 
length and the fluid properties, the Nusselt, Reynolds and Prandtl numbers could be computed from their 
definitions. Meanwhile, the relevant dimensionless parameters for low-speed, forced convection boundary layer 
have been obtained by non-dimensionalizing the differential equation that describes the physical process occurring 
within the boundary layer. An alternative approach is the use of dimensional analysis in the form of BuckingHam 
Pi theorem. The success of the theorem depends on the ability to select from intuition the various parameters that 
influence the problem. Therefore, knowing that 

h = f ( K,C p p,p,V,L) (1.6) 

Before hand one could use the BuckingHam Pi theorem to obtain h, in equation (1.6), [2] 

II. MODELLING AND EXPERIMENTAL VALIDATION 
MODEL FORMATION: To successfully create non-dimensional groups, each time a need arises, a set of rules 
must be followed [3]; the Raleigh method and Buckingham’s Pi-Theorems are reliable. The Raleigh method is an 
elementary technique for finding a functional relationship between variables. Although very simple, the method 
does not provide any information concerning the number of dimensionless group that can be obtained. Another 
drawback of the method is that it can only be used for the determination of the expression for variables that depend 
on a maximum of three or four independent variables only.The BuckingHam Pi theorem is basically an 
improvement over the Raleigh’s method. Apart from its advantage of being able to handle large sets of variables, 
it gives a ready clue on how many dimensionless groups are designated by Pi. 


f(a,h,K,D,Re,P r ) =0 

(2.1) 

„ vVD 


Re - 

(2.2) 

P 


Pr = ^ 

(2.3) 

k 


Therefore 


f{a(h,Kv,D,p,V,p,Cp)} = 0 

(2.4) 


where a. is a constant 

Choosing M, L, T, and K as fundamental dimension for mass, length, time and temperature, implies that: 

Number of Fundamental Dimension, m = 4 

Number of Quantities, n = 7 

Therefore 

Nrj =7-4 = 3 

Since m is 4, there will be four repeating quantities: Geometric (D). Flow (v), fluid density (p) and specific heat 
capacity (Cp) 


Table 2.1: Dimension of Quantities 


S/No 

Quantity 

Symbol 

Unit 

Dimension 

1 

Diameter 

D 

M 

L 

2 

Density 

P 

Kgm" 3 

ML' 3 

3 

Velocity 

V 

m/s 

LT 1 

4 

Viscosity (absolute) 

i 

kgm'V 1 kgms" 

ML-'T' 1 

5 

Thermal Conductivity 

K 

3 K _1 m 2 s 2 K _1 

MLT 3 K-' 

6 

Heat Capacity 

C 

kgs" 3 K' ! 

L 2 T' 2 K 

7 

Heat Transfer Coefficient 

H 


MT- 3 K-' 


K — terms 


n x = D" l V xl p yl Cp zK 

k 2 = D w2 V x2 p y2 Cp z 2 p 
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tt 3 = D w3 V x3 p y3 Cp z 3 h (2.7) 

Setting up and solving equations for n groups, using the dimensions in table 2.1 gives 
For n x = D wX V xX p yX Cp Z X K V 


Putting this dimension terms, gives 

7T X : (LW 1 ) (LT' 1 ) X1 (ML' 3 ) Y1 (L 2 T' 2 K _1 ) = M 0 L°T 0 K 0 

By equating the powers of the dimensions on both sides of equation (2.8), we have 


( 2 . 8 ) 


M: Yi + 1 = O Yi = - 1 

(2.9) 

L: Wi+Xi-3Y, + 2Zi+ 1 = 0 

(2.10) 

T:-Xi - 2Zi -3 = 0 

(2.11) 

K: - Zi -1=0 

(2.12) 

Substituting (2.12) into (2.11) 

-Xi - 2 (-1) -3 = 0 

(2.13) 

-Xi + 2-3 = 0 Xi = - 1 

(2.14) 

Substituting (2.9), (2.12) and (2.14) into (2.10) 

W, - 1 - 3 (-1) +2 (-1)+ 1=0 

(2.15) 

Wi +1-3-2 +1=0 W, = - 1 

Hence, 

(2.16) 

7i x = D wX V xl p yX Cp zK 

(2.17) 

K 

Or 7T X = - 

(2.18) 


pVDCp 


Using similar procedure for equations (2.6) and (2.7), we obtain values for 712 , 713 , as in (2.19) and (2.20) 
respectively. 


Hence 71 2 = 


and 


pVD 

h 


p v Ci 


(2.19) 


( 2 . 20 ) 


Combining equation (2.18), (2.30) and (2.39), we have 


f {a 


kv 




h 


pVDCp vVD pVCp 


= 0 


( 2 . 21 ) 


Hence 


h 


pVCp 


= a 




kv 


pVCp pVDCp 


( 2 . 22 ) 


Using h as the dependent variable and V the independent variable and applying equation 2.18,2.19 and 2.20 to 
obtain the values of K x , 7T 2 and n 3 . 


Therefore, from equation 2.22, 71 2 = a ( 71^ 71 f) (2.23) 

Plotting the independent function ( 71 2 ) against the dependent function ( 7T l ),( 7T 1 ), we obtain graphs which when 
regressed yields the theoretical heat transfer coefficient (hpj. 


\8\ 
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MODEL VALIDATION 

EXPERIMENTAL MATERIALS AND METHOD: Unripe plantains purchased from a local market. Port 
Harcourt, were used in the experiment. The moisture content of the plantain before drying was 26.76g dry base. 
The moisture content was determined manually by periodically weighing of the sample at 3. minutes interval, for 
three hours (61 data points). According to [2], for evaporative heating based on heat and mass transfer analogy, 
as the gas flow over the moist material evaporation occur from the surface, and the energy associated with the 
phase change is the latent heat of vaporization of the liquid. 


Gas 

Flow 


Convective heat 


Evaporative heat 


Fig. Heat and Mass Transfer Analogy 


Applying conservation of energy to a control surface about the material 

d convective d added ~ d evaporation 


(2.24) 


Where q" evaporalio „ may be approximated as the product of the mass (moisture loss) and the latent heat of 


vaporization. 

dx 

d evaporation = nA ^fg = ^ 

(2.25) 

d" convective = hA ( T lf ~ T ,) 

(2.26) 

Since no heat is added, and for constant drying 

d evaporation d evaporation 

Therefore 

(2.27) 

hA (T hf -T) = hA(AA) =^A 

(2.28) 

dx . 



Where — A represents the constant drying rate 
dt 


For surface temperature, T s = 28°C 
and the heating fluid is at T 5 , Tht = T 5 


T S ~T S = AT 5 


Therefore ( h E AAT) 5 = A 


h 

n E( 3 ) 


dx 

dt^ 


Adx 

A(AT)dt j 


When heating fluid is at Ti, Tht = Ti 

\ 

1 


h 

n E( 1 ) 


r Adx ^ 


K A(AT)dt y 

Total Area of Plantain Chips A = 


( nD 1 ^ 




D = diameter of Plantain Chips = 0.033m 
For the 6 pieces of plantain = 0.198m 


(2.29) 

(2.30) 

(2.31) 
(2.32) 

(2.33) 

(2.34) 


191 


Wolume 3 1 Issue 21 


www. ijrtem. com 















Modeling Heat Transfer Coeficient Of Air... 


Hence A 


A ;r0.198 2A 


= 0.03079571 2 


From equation 3.79, 3.80 and 3.81, the mean experimental heat transfer coefficient 


h E — 


A 

A 


dx 

(A T)dt 


\ r 

+ 


J\ 


v 


dx 
(A T)dt 


\ r 

+ 

73 V 


dx 

(A T)dt 


\ 


7 5 


(2.35) 


III. RESULTS AND DISCUSSION 

The results obtained from this work are presented in Table 3.1-3.5 and Figure 3.1-3.5 


Results 


Table 3.1: _ Calibration of the Batch Dryer 


Number of turns of Knob 
(Knob rotation) 

Temperature(°C) 




Velocity(m/s) 




Vi 

v 2 

V A 

V B 

T1 

t 2 

T av 

0 (0°) 

3.60 

3.40 

3.50 

0.66 

66.00 

66.00 

66,00 

3 (1080°) 

4.20 

4.20 

4.20 

0.79 

60.15 

60.05 

60.10 

6(2160°) 

5.00 

4.80 

4.90 

0.92 

54.50 

54.40 

54.45 

9 (3240°) 

5.65 

5.75 

5.70 

1.07 

48.30 

48.20 

48.25 

12 (4320°) 

6.40 

6.40 

6.40 

1.20 

42.00 

42.00 

42.00 


Where V and T are the average values 

Substituting the thermo-physical properties and average velocities based on the calibration, into the 
Dittus-Boelter equation 

hD = — 0023Re “ 8 Pr 04 , (3.1) 

D 

we obtain values of heat transfer coefficient shown in Table 3.1. 


Table 3.2: Thermo-physical Properties and Calculated Values of Heat Transfer Coefficient using Dittus- 

Boelter 


Temperature (°C) 

66.00 

60.10 

54.45 

48.25 

42.00 

Velocity (m/s) 

0.66 

0.79 

0.92 

1.07 

1.20 

Density (kg/m 3 ) 

1.0383 

1.0516 

1.0646 

1.0749 

1.0982 

Specific Heat Capacity (m 2 s- 2 K 1 ) 

1008.3 

1008.2 

1008.0 

1007.8 

1008.5 

Dynamic Viscosity (kgm 's 1 ) 

202. IE 

200.18E-7 

198.4E-7 

196.9E-7 

193.6E-7 

Reynolds Number (-) 

10070.64 

12325.75 

14661.86 

17348.53 

20216.86 

Prandtl Number (-) 

0.702 

0.7022 

0.7028 

0.703 

0.705 

Heat Transfer Coefficient (Wm' 2 K _1 ) 

25.5297 

30.0116 

34.4787 

39.4635 

44.6411 


Modelled Result: The values of the variables of table 3.2 are used to obtain the values of 7T, , 71 2 , 71 ^, as shown 
in Table 3.3. 


Table 3.3: 71 — Groups 


S/No 

1 

2 

3 

4 

5 

7Z y 

0.001121 

0.0009245 

0.0007842 

0.0006681 

0.0005831 


0.0000993 

0.0000811 

0.0000696 

0.0000572 

0.0000495 

^"3 

0.03696 

0.03583 

0.03492 

0.03402 

0.03361 
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Experimental Result: Using the various air flow velocities experiments were performed in order to obtain the 
gradient for drying the unripe plantain. The experimental results have been shown in figures 3.3 to 3.5 where the 
graphs of moisture content versus time at vi, V 3 , and V 5 are shown below. 


n 


0.0375 

0.037 

0.0365 

0.036 

0.0355 

0.035 

0.0345 

0.034 

0.0335 

0.033 


0 


0.0002 


0.0004 


0.0006 

Wj 

Fig 3.1: Graph of u 3 Versus n. 


0.0008 


ir 3 = 6.370x rij + 0.029 
R 2 = 0.997 


0.001 


0.0012 


tt 3 = 68.88HJ + 0.030 

0.0375 -j R 2 = 0.997 



Fig. 3.2: Graph of n3 versus m 


The drying rate (evaporative heat) which is a product of evaporative flux and latent heat of vaporization obtained 
from Figure 3.3, 3.4 and 3.5, and the change in temperature based on Vi, V 2 and V 3 were used to obtain the 
experimental heat transfer coefficient as illustrated below: 

Thf ( i) = 42.00°C, hence A T = 14.00°C 

Thf (3 ) = 54.45°C, hence A T = 26.45°C 

Thf (5) = 66.00°C, hence A T = 38.00°C 

A. = Latent Heat of Vaporization of water = 250 lKJ/Kg [4] 


dx 

dt 


( 1 ) 


dx 

dt 


( 3 ) 


- 0.176g /min = 0.00000285 Kg /s 
-0.173^/min = 0.00000288%/s 


— (5) = -0.17%/min = 0.00000293%/s 
dt 


h 


E 


2501000 

3x0.030795 


" 0.00000285 
v 14.00 


0.00000288 

26.45 


0.00000293^ 
38.00 , 
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h E = 31.6379 Wm~ 2 K l 


Model Heat Transfer Coefficient (Iib) = 34.8524Wm' 2 K _1 
Experimental Heat Transfer Coefficient (Iie) = 31.679Wm" 1 K" 1 

34.8524-31.6379 1/1/w n 

% error = - xl00% =9.22% 

34.8524 

Therefore, a good agreement of about 91% was achieved. 



o *= 


50 
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nT 
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o 
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3 

</> 

O 


Time (Min) 

Fig. 3.3: Moisture Content versus Time (at Vi) 



Figure 3.4: Moisture Content Versus Time (at V 3 ) 


O) 

O 



Time (Min) 

Figure 3.5: Moisture Content Versus Time (at V s ) 

IV. DISCUSSION OF RESULTS 

The values of the velocities and temperatures from the calibration and values from the thermo-physical table 
obtained from literature were substituted into the Dittus-Boelter equation to obtain the relationship between heat 
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transfer coefficients (dependent variables) and velocities (independent variables), which form the basis for the 
Bucking Ham’s Pi-Theorem. 

The BuckingHam’s Pi-Theorem, which uses dimensional analysis, was then used to obtain the Pi groups ( K x , 71 1 

, 77 3 ). Based on the Pi groups obtained, plots of 7T 3 versus 71 x and 71 3 versus 71 2 were obtained and regressed, 

using Microsoft Excel to obtain the modeled equations. The experiment when carried out at velocities of 0.66, 
0.92 and 1.20 ms" 1 , it was observed that the moisture in unripe plantain evaporated faster to approach dryness in 
the order; 0.66ms’, 0.92 ms 1 and 1.2ms' 1 . This is because higher velocities pickup more moisture from the system. 
From the plots of Figures 3.3, 3.4 and 3.5, the equations at the constant drying periods at various velocities were 
also obtained. Finally, a comparison of heat transfer coefficients obtained from the theoretical Buckingham Pi 
Theorem (modeled) and that obtained from the experimental result in the range of velocities investigated showed 
minimal variation of less than 10% 


V. CONCLUSION 

The heat and mass transfer analogy from Newton law of cooling has been shown to be a reliable correlation for 
obtaining heat transfer coefficient experimentally; also proven is that the Buckingham Pi-Theorem is a good and 
simplified method of obtaining correlation from experimental results. The comparison of both the experimental 
and model results shows a percentage error of 9.22%. The variation of the experimental and model heat 
coefficients could be as a result of experimental errors, although it is within acceptable level. 
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